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First-principles study of magnetic properties in Fe-ladder compound BaFe2S3 
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We study the magnetic, structural, and electronic properties of the recently discovered iron- 
based superconductor BaFe 2 S 3 based on density functional theory with the generalized gradient 
approximation. The calculations show that the magnetic alignment in which the spins are coupled 
ferromagnetically along the rung and antiferromagnetically along the leg is the most stable in the 
possible magnetic structure within an Fe-ladder and is further stabilized with the periodicity char¬ 
acterized by the wave vector Q=(7r,7r,0), leading to the experimentally observed magnetic ground 
state. The magnetic exchange interaction between the Fe-ladders creates a tiny energy gap, whose 
size is in excellent agreement with the experiments. Applied pressure suppresses the energy gap and 
leads to an insulator-metal transition. Finally, we also discuss what type of orbitals can play crucial 
roles on the magnetic and insulator-metal transition. 


I. INTRODUCTION 

The discovery of superconductivity in fluorine doped 
LaFeAsO has stimulated intensive studies of iron-based 
superconductors In addition to the high superconduct¬ 
ing transition temperature and the strong flexibility for 
the constituent atoms, similarities of magnetic and struc¬ 
tural properties with the high-Tc cuprate superconduc¬ 
tors are suggestive to understand the pairing mecha- 
nismi^ For instance, the parent material undergoes an¬ 
tiferromagnetic transition under ambient pressure, and 
the superconducting transition appears in the vicinity of 
the antiferromagnetic phase by doping carriers or by ap¬ 
plying pressure. 

Recently, superconductivity in BaFe2S3 has been dis¬ 
covered under pressure with the maximum ^ 14 Ki^ 
Investigation of the magnetic and electronic properties 
in relation to the structural property in this ladder sys¬ 
tem is expected to provide a new path to understand the 
high Tc superconductivity in the iron-based superconduc¬ 
tors. The quasi one-dimensional Fe-ladder formation of 
BaFe2S3, as shown in Fig. [H exhibits a striking differ¬ 
ence from the conventional iron-based superconductors, 
which have two-dimensional Fe networks, but rather re¬ 
minds us of the superconducting two-leg ladder cuprate 
Sri 4 Cu 2404 i<^ While the magnetic structure of BaFe2S3 
has a marked similarity to conventional iron-based su¬ 
perconductors as will be discussed below, BaFe2S3 has a 
semiconducting ground state as contrasted to the metal¬ 
lic ground states of other iron-based superconductors. 

A number of structurally related compounds of 
BaFe2S3, i.e. A¥e2Xz (A=K, Rb, Cs, or Ba and X = S, 
Se, Te), have been discoveredIn particular, the block 
magnetism of BaFe2Se3 has attracted increasing atten¬ 
tion as the novel magnetic structure among iron-based 
materials and has been intensively studied both theoret- 
ically^^— and experimentallyi^d^^— On the other hand, 
the physical properties of BaFe2S3 are reported only re¬ 
cently in association with the discovery of the high-Tc 
superconducting transition^ BaFe2S3 crystallizes in an 


orthorhombic crystal structure whose symmetry belongs 
to space group Cmcm (No. 63 ), which is the same as 
those of CsFe2Se3 and KFe2Se _^^i^^ but is different from 
Pnma (No. 62 ) of BaFe2Se3!^. While the block mag¬ 
netism in BaFe2Se3 is characterized by Q=( 7 r, 7 r, 7 r), the 
periodicity of the stripe magnetism in BaFe2S3 is charac¬ 
terized by Q=( 7 r, 7 r, 0 ). The same stripe-type magnetism 
with BaFe2S3 is also observed in CsFe2Se3 and KFe2Se3. 

The local magnetic moment of iron is reported as 
in BaFe2S3^ which is quite smaller than l.S/is 
of CsFe2Se3^ and 2 . 8 fiB of BaFe2Se3^, and the Neel 
temperatures of 120 K, 175 K, and 255 K for BaFe2S3, 
CsFe2Se3, and BaFe2Se3, respectively, are correspond¬ 
ing to the size of the magnetic moments. The activation 
energy of BaFe2S3 is reported to be less than 0.1 eV^i^ 
which is again smaller than 0.13 eV and 0.178 eV re¬ 
ported for BaFe2Se .j^^i^^ and 0.34 eV for CsFe 2 Se 3 i 2 ^ 

The first-principles study for iron-based superconduc¬ 
tors have provided fruitful information concerning the 
structural, magnetic, and electronic propertiesi^^— For 
BaFe2S3, we recently derived a low energy effective 
Hamiltonian of BaFe2S3 based on a first-principles ap¬ 
proach In this paper, we investigate the magnetic and 
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FIG. 2. (Color online) Possible inequivalent magnetic config¬ 
urations within the Fe-ladder. 


electronic properties in relation to the structural prop¬ 
erty under pressure by using the first-principles calcu¬ 
lations based on the generalized gradient approximation 
(GGA). We show that the stripe spin-ladder magnetic 
order, which is refered to type-I magnetic configuration 
in Fig. [5] with the ordering wave vector Q=(7r,7r,0), is 
the most stable magnetic structure over the investigated 
range of pressure. We also discuss how the electronic 
structure changes under pressure and undergoes a metal- 
insulator transition, which is necessary for the supercon¬ 
ducting transition. 


ting the 2x2x1 super cell, corresponding to the magnetic 
unit cell for the periodicity characterized by Q=(7r,7r,0). 
The A:-point mesh is set as the (2,2,2) for the self- 
consistent calculations and (8,8,8) for plotting the den¬ 
sity of states (DOS) for the calculations of the super 
cell. The calculations were performed with the experi¬ 
mental lattice constants a=8.78 A, 6=11.23 A, c=5.29 A 
for ambient pressure, and the pressure effect are consid¬ 
ered as change in the lattice constants observed exper- 

e=(0,0,0) 



FIG. 3. (Golor online) Gomputed spin distribution of the 
type-I magnetic order characterized by Q=(0,0,0), (7r,0,0), 
(tTjTTjO). The yellow (blue) surface is an isosurface of positive 
(negative) spin density. 


II. METHOD 

The first-principles calculations are implemented in 
the VASP program code with projector augmented-wave 
(PAW) method using the exchange-correlation functional 
proposed by Perdew, Burke, and Ernzerhofi^ We sys¬ 
tematically explored all the possible spin configurations 
within the Fe-ladder, as shown in Fig. [2l for the char¬ 
acteristic wave vectors Q=(0,0,0), (7r,0,0), (7r,7r,0), set- 


imentallyi^ The atomic relaxation is performed for the 
sulfur atoms within the x-y plane, starting from the ini¬ 
tial atomic configuration as the experimental values i.e. 
Ba(4c) (0.0,0.686,0.25), Fe(8e) (0.154,0.0,0.0), S(4c) (0.0, 
0.II6, 0.25), and S(8g) (0.208, 0.378, 0.25), using the ex¬ 
perimental lattice constantsi^ The energy cutoff for the 
plane waves is set as 500 eV through the calculations. 
The convergence conditions for the total energies are set 
as 10“^ eV for electronic self-consistent loop and 10“^ 
eV for the ionic relaxation. 


III. RESULTS 


netisin with different Q are shown in Fig. |31 Galculated 


Firstly, we discuss the stable magnetic structure of 
BaFe 2 S 3 . The magnetic distributions of type-I mag- 
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Experimental structure (P=0) Optimized structure (P=0) Optimized structure (P=16.27 GPa) 

Total energy (eV) Fe-moment (/is) Total Energy (eV) Fe-moment (/is) Total energy (eV) Fe-moment (/is) 


Q=(0,0,0) 

Q=(^,0,0) 

Q=(7r,7r,0) 

-291.006 

-290.786 

-291.050 

1.28 

1.08 

1.03 

-293.901 

-293.901 

-293.901 

0.00 

0.00 

0.00 

-293.296 

-293.295 

-293.296 

0.00 

0.00 

0.00 




Type-I 




Experimental structure (P=0) 

Optimized structure (P=0) 

Optimized structure (P=16.27 GPa) 

Total energy (eV) Fe-moment (/is) Total Energy (eV) Fe-moment (/is) Total energy (eV) 

Fe-moment (/is) 

Q=(0,0,0) 

-296.009 

2.26 

-296.778 

2.08 

-293.558 

1.42 

Q=(^,0,0) 

-295.967 

2.27 

-296.742 

2.08 

-293.632 

1.10 

Q=(7r,7r,0) 

-296.045 

2.27 

-296.809 

2.07 

-293.672 

1.40 




Typ( 

5-II 




Experimental structure (P=0) 

Optimized structure (P=0) 

Optimized structure (P=16.27 GPa) 

Total energy (eV) Fe-moment (/is) Total energy (eV) Fe-moment (/is) Total energy (eV) 

Fe-moment (/is) 

Q=(0,0,0) 

-294.964 

2.55 

-295.245 

2.42 

-293.296 

0.00 

Q=(^,0,0) 

-294.953 

2.49 

-295.310 

2.31 

-293.296 

0.00 

Q=(7r,7r,0) 

-294.653 

2.49 

-295.043 

2.30 

-293.296 

0.00 




Type-III 




Experimental structure (P=0) 

Optimized structure (P=0) 

Optimized structure (P=16.27 GPa) 

Total energy (eV) Fe-moment (/is) Total Energy (eV) Fe-moment (/is) Total energy (eV) 

Fe-moment (/is) 

Q=(0,0,0) 

-294.402 

2.15 

-295.314 

1.87 

-293.295 

0.15 

Q=(^,0,0) 

-294.259 

2.14 

-295.182 

1.89 

-293.296 

0.00 

Q=(7r,7r,0) 

-294.281 

2.15 

-295.195 

1.86 

-293.295 

0.14 


TABLE I. Total energies and spin moments calculated for the different spin configurations (see text). 


total energies and Fe magnetic moments are shown in ta¬ 
ble ID As shown in the table, the calculations predict the 
most stable magnetic structure as the type-I magnetic 
alignment, in which the spins are coupled ferromagneti- 
cally along the rung and antiferromagnetically along the 
leg, with the periodicity characterized by the Q=(7r,7r,0), 
which corresponds to the experimentally observed mag¬ 
netic alignment. The energy difference is mainly dom¬ 
inated by the spin configuration within the Fe-ladders, 
changing the total energy with the order of 1 eV, and 
the different Q structure change it with the order of 0.1 
eV. 

Relaxation of sulfur atoms somewhat reduces the mag¬ 
netic moments from those for the experimental atomic 
position for the type-I, II, and III magnetic states and 
completely suppresses it for FM spin alignment. Remark¬ 
ably, the type-I magnetic configuration is the most stable 
magnetic structure within an Fe-ladder in spite of the 
smaller magnetic moments compared to those of type- 
II. This result indicates that the atomic configuration 
strongly affects to stabilize the type-I magnetism through 
coupling with the magnetic moments. The predicted Fe’s 
magnetic moments 2.07^b of the type-I magnetic order 
is considerably larger than the experimental value I.S/tb, 
and this kind of overestimation for Fe’s moments was also 
reported for LaFeAsO. 36,37 

We also investigated the magnetic states under pres¬ 


sure by using the experimental lattice constants with per¬ 
forming relaxation of sulfur atoms. As shown in table HI 
whereas the magnetic moments of the type-I magnetic 
states are not much reduced at P=16.27 GPa, those of 
the other magnetic states are strongly suppressed under 
the high pressure. To examine the pressure dependence 
of the magnetic states in further detail, we show in Fig.|4] 
the detailed pressure dependence of the total energies and 
Fe-spin moments for the magnetic states with the peri¬ 
odicity Q=(7r,7r,0). From the figure we can find the clear 
difference in the pressure dependence for all of the mag¬ 
netic structures, though the type-I, II, and III magnetism 
exhibit the same order of the magnetic moments at ambi¬ 
ent pressure. The pressure effect strongly suppresses the 
magnetic moments of the type-II magnetism, in which 
the magnetic moment vanishes already at P=2 GPa and 
gradually suppress the magnetic moments of the type-III 
magnetism, leading to the nonmagnetic state around the 
16 GPa. On the other hand, the pressure dependence of 
the magnetic moment for the type-I magnetism is rather 
weak and leave the magnetic order stable even under high 
pressure. 

The optimized sulfur positions, whose coordinates are 
defined in Fig. [SJ are ?/si=l-18 A, a;s2=a;s2'=2.52 A, and 
yS2= 2/S2'=1-36 a while the experimental positions are 
?/si=l-30 A, a;s2=a;s2'=2.57 A, and ?/S2=ys2'=l-37 A. 
The deviation between the calculation and the experi- 
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FIG. 4. (Color online) Pressure dependence of (a) total energy 
and (b) magnetic moments for all possible magnetic structures 
characterized by Q = (7r,7r,0). 


merit for the SI position is somewhat larger than that for 
the S2 atom, yet the deviations are about 1% at most for 
the lattice constants. We also show the pressure depen¬ 
dence of the sulfur’s atomic positions in Fig. 0 in which 
the positions of sulfurs are measured from the Fe-ladder 
as defined in Fig. [S] (a). Interestingly, from Fig. [S] (b), 
significant shifts in the sulfur positions are only observed 
for the sulfurs out of the Fe-ladders (S2 and S2’), and the 
positions of sulfurs located between the ladder’s legs (SI) 
does not strongly depend on the pressure. Note that S2 
and S2’ in Fig.^are equivalent in the non-magnetic crys¬ 
tal structure but are inequivalent in the magnetic state. 
The calculations indeed show small deviation of the two 
sulfur positions above P ~12 GPa as shown in Fig.[S](b). 

Finally, we discuss the electronic structure of BaFe 2 S 3 . 
Experimentally, BaFe 2 S 3 is semiconductor with the small 
energy gap of 0.06-0.07 eV As shown in Fig. |6l the 
energy gap of the BaFe 2 S 3 is not produced with the 
Q=(0,0,0) periodicity but produced with the periodicity 
of Q=(7r,7r,0) with the size 0.07 eV for the type-1 mag¬ 
netism, which is in excellent agreement with the experi¬ 
mental evaluations. This result indicates that the mag¬ 
netic exchange coupling between the Fe ladders is crucial 
to create the tiny energy gap. Figure [7] shows pressure 
dependence of total DOS for the Q=(7r,7r,0) type-1 anti¬ 
ferromagnetic states around the Fermi level. The pres- 
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FIG. 5. (Color online) (a)Definitions of sulfur’s relative po¬ 
sitions and (b) the calculated values for the type-I magnetic 
states. The red surface sheet in figure (a) indicates a surface 
on which the Fe-ladder lies. 


sure strongly suppresses the energy gap and leads to a 
metallic state around P ~5 GPa. The peak structure 
of DOS located below the Fermi level comes close to the 
Fermi level as the pressure increases and passes over the 
Fermi level at P ~12 GPa. As a result, pressure strongly 
enhances the DOS at the Fermi level which takes its max¬ 
imum at P ^12 GPa. The orbital resolved DOS of Fe-d 
orbitals for the pressure P=0 and 12.36 GPa is shown in 
Fig. [H The large difference in the structure of the DOS 
for the majority and minority spin indicates that all of 
the Fe-d orbitals contribute to the magnetic order. More 
precisely, contribution for the spin-moment of each d or¬ 
bital calculated from the spin-DOS is 0.37 (0.28) for d^y, 
0.39 (0.30) for dy^, 0.34 (0.27) for d 3 ^ 2 _^ 2 , 0.47 (0.36) 
for dxz, and 0.49 (0.37) for leading to the to¬ 

tal Fe-d moment 2.06 (1.58), at P=0 (12.36 GPa). The 
larger magnetic moments of the dxz and dx^-y'^ orbitals 
compared to other d-orbitals show the primary contri¬ 
bution of these two orbitals for the magnetism, which is 
consistent to our recent studies showing that low-energy 
effective Hamiltonians can be constructed from these two 
orbitalsi^ The d^x'^_x 2 orbital mainly contributes to the 
DOS peak below the Fermi level at ambient pressure but 
applied pressure suppresses it as in Fig.[8](b) and (d). Al¬ 
ternatively, the dyz and the dx^-y^ orbitals develop with 
applied pressure and contribute to the maximum DOS 
peak at Fermi level at P ~12 GPa. Although further in¬ 
vestigation is necessary to reveal the orbital contribution 
to the superconducting transition, the interplay of these 
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FIG. 6. (Color online) Density of states of Type-I AFM 
magnetic ordered states with the periodicity characterized by 
Q=(0,0,0) and by Q=(7r,7r,0). Insets are the enlarged views 
around the Fermi levels. 



Energy (eV) 

FIG. 7. (Golor online) Pressure dependence of density of 
states for type-I AFM magnetic ordered states with the peri¬ 
odicity of Q=(7r,7r,0). 


temperature. 


IV. CONCLUSION 

We investigated the magnetic, structural, and elec¬ 
tronic properties of BaFe 2 S 3 based on the first-principles 
calculations. The GGA calculations predicted the stripe- 
type spin-ladder magnetism, whose magnetic periodic¬ 
ity is characterized by Q=(7r,7r,0), in the ground state. 
The calculations of the pressure dependence for the to¬ 
tal energy and magnetic moments indicate that the most 
stable type-I magnetic structure is not sensitive against 
the applied pressure compared to other magnetic struc¬ 
ture and leave the magnetic order stable even at high 
pressure as a result. The magnetic exchange interac¬ 
tion between the Fe ladders creates a tiny energy gap 
at ambient pressure and is therefore a crucial ingredi¬ 
ent to reproduce the insulating magnetic ground state as 
observed experimentally. Applied pressure strongly sup¬ 
presses the energy gap and leads to a metallic magnetic 
state, bringing the high DOS at the Fermi level due to the 
formation of the DOS peak structure around the Fermi 
level. The magnetic moment predicted at the ambient 
pressure is considerably larger than that observed exper¬ 
imentally. The discrepancy of the magnetic states may 
be improved with more advanced first-principles method 
considering the strong electron correlation in the iron- 
based compound and more precise prediction of the 
magnetic states is our future issue. 



Fe-d orbitals should be crucial for the high-Tc transition 


FIG. 8. (Color online) Density of states of Fe-d orbitals for 
Type-I (Q=(7r,7r,0)) AFM states for P—0 and 12.36 GPa. 
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